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We consider diagonal matrix elements of local operators between multi-soliton states in finite 
volume in the sine-Gordon model, and formulate a conjecture regarding their finite size dependence 
which is valid up to corrections exponential in the volume. This conjecture extends the results of 
Pozsgay and Takacs which were only valid for diagonal scattering. In order to test the conjecture we 
implement a numerical renormalization group improved truncated conformal space approach. The 
numerical comparisons confirm the conjecture, which is expected to be valid for general integrable 
field theories. The conjectured formula can be used to evaluate finite temperature one-point and 
two-point functions using recently developed methods. 
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1. INTRODUCTION 

In this paper we continue a program to describe form 
factors in finite volume, initiated in [1, 2] and further 
developed and extended in the works [3-7]. Previously, 
finite volume form factors have also been studied in other 
approaches [8-10]. 

The present effort at describing finite volume form fac- 
tors is directed towards an understanding of finite tem- 
perature correlation functions and other related cases 
where the spectral expansion of a quantity is ill-defined 
due to singularities associated with disconnected contri- 
butions present in operator matrix elements. A different 
attempt to find a spectral expansion for finite tempera- 
ture correlators is Doyon's finite temperature form factor 
formalism [11]. 

We work in two-dimensional integrable quantum field 
theories, where the S matrix can be obtained from the 
bootstrap [12]. In such theories, using the scattering am- 
plitudes as input it is possible to obtain a set of equations 
satisfied by the form factors [13]. The complete system 
of form factor equations, which provides the basis for the 
form factor bootstrap, was proposed in [14]. For a de- 
tailed and thorough exposition of the subject we refer to 
[15]. 

Finite volume form factors can be used as a tool to 
compute finite temperature correlation functions [2, 16- 
18], or matrix elements of local fields in a boundary set- 
ting [19]. Another application of the formalism is the 
extension of form factor perturbation theory, which was 
developed in [20] to described non-integrable perturba- 
tions of integrable quantum field theories, to higher or- 
ders in [21]. 

One of the applications of our framework is that it al- 
lows for a direct comparison of solutions of the form fac- 
tor axioms to field theory dynamics. In fact, this works 
the other way as well, since we can use well-established 



results about the form factor solutions in a given theory 
to test our ideas or conjectures that arise in developing 
the finite volume form factor program. There is still place 
for developments, as we do not yet have a complete de- 
scription of finite volume form factors for the case of non- 
diagonal scattering. The first steps were taken in [5] with 
a study of sine- Gordon breather and two-soliton form 
factors; later it was extended to multi-soliton states [7], 
using the framework for numerical evaluation of multi- 
soliton form factors developed in [22]. Another important 
direction is to incorporate exponential finite size correc- 
tions: so far this was only done for the so-called ^-terms 
[4, 6] , and even in that case the method is not yet entirely 
systematic. 

Theories with non-diagonal scattering, in which the 
spectrum contains some nontrivial particle multiplets 
(typically organized into representations of some group 
symmetry), such as sine-Gordon or the 0(3) nonlinear 
sigma model are very important for condensed matter 
applications (e.g. to spin chains or one-dimensional elec- 
tron systems; for a review see [23]). The finite volume 
description of form factors can be used to develop a 
low-temperature and large-distance expansion for finite- 
temperature correlation functions, which could in turn 
be used to explain experimental data, e.g. from inelastic 
neutron scattering [16, 24]. Another interesting applica- 
tion is to extend the computation of one-point functions 
of bulk operators on a strip to the non-diagonal case. Fi- 
nite volume methods are also a promising approach in 
the study of quantum quenches in integrable quantum 
field theories [19, 25]. 

In this paper we treat the sine-Gordon model as an 
example. It can be considered as the prototype of a 
non-diagonal scattering theory, and it has the advantage 
that its finite volume spectra and form factors can be 
studied very effectively numerically using the truncated 
conformal space approach (TCSA), originally developed 
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by Yurov and Zamolodchikov for the scaling Lee- Yang 
model [26], but later extended to the sine-Gordon the- 
ory [27]. Its exact form factors are also known in full 
generality [15, 28-31], and so it is a useful playground to 
test our theoretical ideas on finite volume form factors. 
However, before embarking on the present program, a 
technical problem had to be solved. Namely, it was clear 
from our earlier studies [5, 7] that truncated conformal 
space did not converge very well for diagonal matrix ele- 
ments. In order to solve this problem we implemented the 
numerical renormalization group (NRG) improvement in- 
troduced by Konik and Adamov [32] . As we demonstrate, 
the resulting NRG-TCSA method proved to be accurate 
enough to perform a stringent test of our conjectures. 

The paper is organized as follows. After a brief review 
of the necessary facts about sine-Gordon model and its 
finite volume soliton form factors in Section 2, we formu- 
late our main conjecture in Section 3. In Section 4 we 
give a brief description of the numerical methods, and 
then present the results of our computations in Section 
5. Section 6 is reserved for the conclusions. 



2. SOLITONS IN FINITE VOLUME 
2.1. Action and S matrix 

Sine-Gordon model is defined by the classical action 

A = J d 2 x Qd^d^ + cos /3$ 

The spectrum of the quantum theory is generated by a 
doublet of a soliton and an antisoliton, both of mass M. 
Their exact 5 matrix can be written as [12] 

s££(o,t) = si%(.e,t)s (e,t) (2.1) 

? 8^ - f3 2 
where the non-zero elements are 



and 



S T (0,S) 



5||(6 


,0 = 


5+1(6 


,0 = 


5++(6 


,0 = 


sinh 


(f) 






sinh ( 


iTT-e\ 
« ) 



i sin 


(1) 




sinh 


{<¥) 



f r°° tit sinh Zil—£>1 
y J \ Jo t sinh 2|* cosh f 



sin 6t 



Besides the solitons, the spectrum of theory contains also 
breathers which are bound states of a soliton and anti- 
soliton. 



Sine-Gordon model can also be represented as a free 
massless boson conformal field theory (CFT) perturbed 
by a relevant operator, with the Hamiltonian 

H = J dx^: (<9 t $) 2 +(d x $) 2 : +A J dx : cos/3$ : (2.2) 

where the semicolon denotes normal ordering in terms of 
the modes of the A = massless field. Due to anoma- 
lous dimension of the normal ordered cosine operator, the 
coupling constant [i has dimension 



A ~ [mass] 



2-/3 2 /4tt 



so it defines the mass scale of the model and the dimen- 
sionless coupling parameter is (3. 

2.2. Soliton form factors 

The class of operators we consider consists of exponen- 
tials of the bosonic field Their vacuum expectation 
value is known exactly [33]: 



x exp 




a 2 f3 2 _ 2 
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with M denoting the soliton mass related to the coupling 
A defined in via [34] 



2T(A) / v^ r (t^a 
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Multi-soliton form factors, i.e. 



F 



H...l N V."l) ■ 



,,0 N ) = (Q\O\A lN (6 N )...A n (0i)), (2.5) 



(A± denoting a soliton/antisoliton and the ordering of 
rapidities 6>i > ... > On) have been constructed using 
several different approaches: the earliest construction is 
by Smirnov (reviewed in [15]), then a free field repre- 
sentation by Lukyanov [28, 29], and later in the work 
by Babujian et al. [30, 31]. Here we use formulae from 
Lukyanov's work [29] (in conjunction with its numeri- 
cal evaluation method given in [22]); however, certain of 
his conventions are different and therefore we change the 
labeling of the form factors accordingly (see eqn. (2.10) 
below) . The reason is that the form factors we use satisfy 
form factor bootstrap relations which are slightly differ- 
ent from Lukyanov's conventions; in this we conform to 
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the conventions of the papers [1, 2]. In our notations, the 
form factor equations are: 

I. Lorentz-invariance 

Ft. iN (0i + a, • • ■ , e N + A) = c s (°) A ^... ljv (0 1 , ...,e N ) 

(2.6) 

where s(O) is the Lorentz spin of the operator O, 

II. Exchange: 

F i 1 ...i k i k+1 ...i N (8ii • • • > 6k, 0fc+i) • • • , 0jv) = 

xi ^?...j fc+ ij fc ...iiv • • ■ > k+i,Ok, ■■■,9n) (2-7) 

III. Cyclic permutation: 

F° i2 .. AN (6 1 +2in,6 2 ,...,6 N ) = 

e 2vMO) Fg_ ANil (e 2 ,...,e N ,9 1 ) (2.8) 

where w(0) is the mutual locality index between the op- 
erator O and the asymptotic field that creates the soli- 
tons. 

IV. Kinematical singularity 

-iRes i FP kil __ AN (e + i7r,e',e 1 ,...,e n ) = 

8=0 

CW - e 2 ™Ws k k ^(9> e 1 )S k k f 2 {9 1 2 ) 

■■■^ in {e'-e N ))Fl„ jN {e u ...,e N ) (2.9) 

where C is the charge conjugation matrix. 

There is a further equation that relates form factors 
containing breathers to those containing only solitons, 
but it is not needed for multi-soliton states. These equa- 
tions are supplemented by the assumption of maximum 
analyticity (i.e. that the form factors are meromorphic 
functions which only have the singularities prescribed by 
the axioms) and possible further conditions expressing 
properties of the particular operator whose form factors 
are sought. 

The form factors of the operator 

O a = e ia ^ 

which satisfy equations (2.6-2.9) with the locality index 

ui(O a ) = a mod 1 
can be obtained from 

F 2i...*m (01) • • • > #2n) = (-l) n ^'_ 2 2JV ..._ <T i (027i! • • • , #l) 

= (-i)"4;;.U(02n,...,0i) 

(2.10) 

where the functions T (originally derived by Lukyanov 
in [29]) are specified in appendix A of [7]. 



2.3. Soliton form factors in finite volume 

The formulae for finite volume form factors, derived 
in [1, 2], were generalized for the case of non-diagonal 
theories in [5] and further investigated in [7]. Here we 
only recall the necessary facts; for more details the reader 
is referred to the original papers. 

In finite volume L, the space of multi-soliton states can 
be labeled by momentum quantum numbers ,Iff. 
We introduce the following notation for them: 

\{h,I 2 ,...M)P (2.11) 

where the index r enumerates the eigenvectors of the n- 
soliton transfer matrix, which can be written as 

T{d\{e 1 ,...,e N \tX = 

(a - ejs^i (v-6 2 )...s:i» iiN (#- 6 N ) 

where 6±, . . . , Off are particle rapidities. The transfer ma- 
trix can be diagonalized simultaneously for all values of 

rm{e 1 ,...,e N })ii:X = 

*w (#,{0k})M:i in ae k }) 

We can assume that the wave function amplitudes ^( r > 
are normalized and form a complete basis: 

E <U *t.i N (mr = s rs (2.i2) 

E *£U *ili.v « W = ■ ■ • <w 

r 

these eigenfunctions describe the possible polarizations of 
the N particle state with rapidities 0\ , . . . , Off inside the 
2 N dimensional internal space indexed by i\ . . . iff. The 
transfer matrix can be diagonalized using the algebraic 
Bethe Ansatz (cf. Appendix A of [5]), which enables 
one to compute the exact form of eigenvalues and 
eigenvectors \Jf( r ). 

The rapidities of the particles in the state (2.11) can 
be determined by solving the quantization conditions 

Qj^(0i, ■ • ■ , 0n) = 

MLsinh0j+^ r) (0i,...,0iv) = 2%I h k = l,...,N 
sl r) (6 1 ,...,6 N ) = -ilog&H6 j ,{6 k }) (2.13) 

When considering rapidities which solve these equations 
with given quantum numbers 1\ , . . . Iff and a specific po- 
larization state r, they will be written with a tilde as 
0i, ... , Off . 

Using the above ingredients, the finite volume matrix 
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elements can then be written as [5, 7] 

^({i' 1 ,...,i' M }\o(0M{h,...,iN}) { [ ) 



0{e~^ L ) (2.14) 



where p( r ) and denote the density of states of types r 
and s, which can be calculated as the Jacobi determinant 
of the Bethe-Yang equations (2.13), considering them as 
a mapping from the rapidity to quantum number space: 



Furthermore, 



{dQ\ r) 

•' M=dct n^r 



(2.15) 



j,k=l,...,N 



(6' M , . . .,9[\6i, 



'N 



E E *£U mr 



31---3M ll—tN 

xF-° , . . (6' M + in, . . . , 6\ +i7r,0i, 



(r) 



({**}) 



(2.16) 



is the (s, r)-polarized form factor (the bar denotes the 
antiparticle) . The absolute value in (2.14) is necessary to 
account for the different phase conventions of the multi- 
particle states used in the form factor bootstrap and in 
the finite volume calculations. 



3. A CONJECTURE FOR DIAGONAL MATRIX 
ELEMENTS 

Relation (2.14) is only valid for matrix elements with 
no disconnected pieces, i.e. when the rapidities in the 
two finite volume states are all different from each other. 
If there are particles with exactly coinciding rapidities 
in the two states, i.e. 0' k = 0i for some k and I, then 
there are further contributions. Note that equality of 
two quantum numbers such as I' k = Ii is not sufficient 
for the presence of a disconnected contribution, as the 
corresponding rapidities will in general be different due 
to the terms involving the phase shifts 6j . Therefore 
such terms are only present for the case when the two sets 
of quantum numbers are exactly identical, and also in the 
special case when the two states each contain a particle 
with exactly zero rapidity. At present, the disconnected 
terms are only known for states with diagonal scattering; 
the form of these contributions was obtained in [2] . 

Here we present a conjecture for diagonal matrix el- 
ements with non-diagonal scattering, which is an edu- 
cated guess based on the results valid for diagonal scat- 
tering and also on some lessons learned from nested Bethe 
Ansatz systems [35, 36]. For diagonal scattering, the for- 



mula introduced in [2] states that the diagonal finite vol- 
ume matrix element can be computed as 

({h,...,I N }\O(0,0)\{h,...,I N }) L = 

1 



]T F(A) L p{A) L +0(e-» L ) 



Ac{l,2,...JV} 



(3.1) 



where A = {1,2,... N} \ A. Let us denote the elements 
of the sets A and A as 

A = {A 1 ,...,A l ] 
A = {A u ...,A N -i} 

with I = \A\ the cardinal number (number of elements) 
of the set A. Then 

p{A) L = p{6 Al ,...,6 AN _ l )L 

is the Z-particle Bethe-Yang Jacobi determinant involving 
only the subset A of the N particles, and 



r(A) L = Ff(e Al ,...,e Al ) 



where 



Ff(0i, . . .,61)^. 



lim F U {6 1 + in + e, . . . , 6 1 + m + e, U . . . , l ) f fli 



is the so-called symmetric evaluation of the diagonal form 
factor involving the particles in set A (the bar denotes 
the antiparticle). 

The main observation is that the density p is only well- 
defined for states that have an internal state which is an 
eigenvector of the multi-soliton transfer matrix. There- 
fore, in order to define the disconnected term correspond- 
ing to a given subset A, the 7V-soliton wave function am- 
plitude must be decomposed accordingly. Let us suppose 
that 

\{I 1 ,I it ...,I N })% ) 

is a finite volume state which corresponds to a wave func- 
tion eigenvector 



({Ok}) 



Using the orthogonality and completeness of the ^ ampli- 
tudes (eqn. (2.12)), for any subset Ac {1, 2, ... N} we 
can define appropriate branching coefficients C to split 
the wave function into pieces given by tensor products 
of transfer matrix eigenvectors for the sets A and A as 
follows: 



Kli* ({**}) = E C # (WW 



5 



where the sum over s runs over all possible polarization 
states (transfer matrix eigenvectors) for I particles, and 
the sum over t runs similarly for polarization states of 
N — I particles. From (2.12) it follows that the branching 
coefficients are normalized as 



s.t 



1 



The conjectured generalization of (3.1) is then 
V({h,...,I N }\O(0,0)\{h,...,Iir})£ ) = 



p(r)({l,...,N}) L 

£|<#({fc} 14 



where 

F {S \A) L = limF°( s 



E 

Ac{l,2,...N] 
2 



(3.2) 

^ s \A) L p^(A) L + 0(e-t iL ) 



\(s) 



is the symmetric diagonal limit of the (s, s)-polarized 
form factor from (2.16) and 

P ®(A) L = pM(e Al ,...,0 Ajf _ l ) 

is the density of states in the channel. 

It is easy to see that eqn. (3.2) reproduces (3.1) for 
diagonal scattering, as in that case there is only a single 
polarization state for any given number of particles and 
the branching coefficients C are all equal to 1. Indeed 
the structure of (3.2) is very obvious, the only ambiguity 
is which evaluation of the diagonal form factor to use; 
for that we substituted the symmetric one, to keep the 
correspondence with (3.1). 

The first nontrivial case arises when N = 2. In this 
subspace there are two states with non-diagonal scatter- 
ing, which contain a soliton and an antisoliton. The wave 
function amplitudes of the transfer matrix eigenstates are 



*(+) = 



s/2 



( ° \ 
1 



which are nothing else but the branching coefficients of 
these states in a (one-soliton) x (one-soliton) basis; it is 
a great simplifying feature that these are rapidity inde- 
pendent. The general formula (3.2) specified to this case 
yields 



W({I 1 ,J 2 }|O(0,0)|{J ll I 2 })i ±) 



,(±) 



limif (±) (0 2 + e A + e\9 1 ,6 a ) 



(±)_ 



(pi(0i)+pi(0 2 )) +p {±) (6 1 ,8 2 )(0) 



(3.3) 



where 



pi(0) = ML cosh 6 



is the one-soliton state density, the diagonal one-soliton 
form factor is 



17T, I 



F 



«7T, I 



and it is independent of 9 due to Lorentz invariance. 
In addition, the (±, ±)-polarized two-soliton-two-soliton 
form factors are given by 



F° 



{6' 2 + m,0' 1 +i<7r,0i,6 2 ) 



rF^ 
sF< 



ITT 



ITT, 



+ - \P2 + * 7r i "1 + '' 7r 

rsF°A 



h) 



i 



with r, s = ±1. The rapidities can be obtained by 
solving the quantization conditions: 

Qf L \0i,6 2 )= ML sinh 9\ + 5±{6 l - 6 2 ) = lixh 
Q 2 ±] (Si^2) = MLsmhd 2 + 5 ± (9 2 - 9 X ) = 2nl 2 (3.4) 

where the phase-shifts S± are defined from the eigenval- 
ues of the two-soliton 5-matrix in the neutral subspace 
by 



S_(9)=S+l(9)-S+l(9) = e 

where the minus sign introduced in the first line is a redef- 
inition, which ensures that the phase-shifts are odd and 
continuous functions of the rapidity 9. Due to this con- 
vention the even states are quantized with half-integer, 
while the odd states are quantized with integer quantum 
numbers. The densities of neutral two-soliton states 
can be written as the Jacobi determinant [5] 



9Q\ ' 9Q\ 



(±) 



,(1) 



902 



902 



4. NUMERICAL METHODS 

4.1. TCSA in a nutshell 

To evaluate the form factors numerically, we use the 
truncated conformal space approach (TCSA) pioneered 
by Yurov and Zamolodchikov [26] . The basic idea of this 
approach is to consider the field theoretic Hamiltonian 
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on the space of states of the conformal field theory, and 
truncate the basis of this space to finitely many vectors 
by placing an appropriate upper energy cutoff. The ex- 
tension to the sine-Gordon model was developed in [27] 
and has found numerous applications since then. The 
Hilbert space can be split by the eigenvalues of the topo- 
logical charge Q (or winding number) and the spatial 
momentum P, where the eigenvalues of the latter are of 
the form 

2ns 

17 

s is called the 'conformal spin'. In sectors with vanishing 
topological charge, we can reduce the size of the Hilbert 
space using the symmetry of the Hamiltonian under con- 
jugation of the solitonic charge: 

To fix our notations, we give here the formal definition 
for the cut-off Hilbert space: 

^s,g ( e c«t) = 

{|*> : Q|*> = g|tt>, (L - Lo)|tf) = s|*>, C|tf) = ±|tf) 
and (L + L - 1/12) I*) = e|*) with e < e cut ) 

where |^) runs over the Hilbert space of the ultraviolet 
massless boson with quasi-periodic boundary conditions 

2ir 

$(x + L,t) = $(x,t) + — q qeZ 

P 

Lq and Lq are the usual Virasoro generators, and e cu t is 
the (dimensionless) cut-off parameter. The upper index 
± corresponding to projections under C is applicable only 
in the Q = sectors. 

For the purpose of form factor calculations, we choose 
the operator 

O =: e i/3$ : 

where the semicolons denote normal ordering with re- 
spect to the A = free massless boson modes. It has the 
conformal dimensions 

B 2 

A = A = £- 
Sir 

Using relation (2.4) we can express all energy levels and 
matrix elements in units of (appropriate powers of) the 
soliton mass M, and we also introduce the dimensionless 
volume variable I = ML. The general procedure is the 
same as in [1, 2]: the particle content of energy levels 
can be identified by matching the numerical TCSA spec- 
trum against the predictions of the Bethe-Yang equations 

(2.13) . After identification, one can compare the appro- 
priate matrix elements to the theoretical values given by 

(2.14) and (3.2). There are some technical issues in the 
identification of states due to level crossings, which can 
also affect numerical accuracy; for a discussion of these 
we refer to [3, 5]. 



4.2. The numerical renormalization group in TCSA 

As mentioned in the introduction, our earlier studies 
showed [5, 7] that TCSA did not converge very well for di- 
agonal matrix elements. This means that the numerically 
evaluated matrix elements showed a marked dependence 
on the cutoff e cut . It is also known from earlier expe- 
rience that TCSA converges faster for smaller values of 
/3 (or £), and this was also clear from previous calcula- 
tions [5, 7, 21, 37]: decreasing £ the cut-off dependence 
was reduced and the agreement between the numerical 
TCSA results and the predictions of the finite volume 
form factor formalism was improved at the same time. 

Following the proposal of Konik and Adamov [32] , we 
use a Wilsonian type numerical renormalization group 
(NRG) to improve the precision of TCSA. This allows 
us to use much higher values for the cutoff than for the 
usual TCSA. Usually, the attainable cutoff is in the range 
e cut = 20 ... 26 (depending on the sector and the value of 
0) , which corresponds to an upper limit of around 20000 
states. E.g. in the Q = 2, s = sector at £ = |, e cut = 26 
translates to 23771 states. With the NRG improvement, 
it is possible to take into account the effect of several 
hundred thousand states, e.g. in the previous example 
we could take into account 840000 states corresponding 
to the cutoff value of 38 (the number of states increases 
exponentially with the cutoff). 

Our implementation of the NRG in TCSA is as follows. 
Let us take a given value e cu t of the cutoff with N states, 
which we order by increasing value of their conformal 
energy eigenvalue. We start with the lowest lying Nq of 
them (which is chosen to correspond to some lower value 
e cu t,o °f the cutoff ). We split the remaining states into 
k "shells", each containing Ni states and a "remainder 
shell" of N r states such that 

N = N + kN 1 + N r 

The idea is that first we diagonalize the lowest Nq x No 
block of the Hamiltonian and retain the first n lowest 
lying eigenvalues with their corresponding eigenvectors. 
Then we add the first shell to the basis and recompute 
the lowest n eigenvalues and corresponding eigenvectors, 
which are now Nq + Ni dimensional. In the next step, 
we add to these modified eigenvectors the next shell, and 
recompute the eigenvalues and eigenvectors, which in the 
original conformal basis will now have Nq + 2N± dimen- 
sions. Repeating the procedure k times and applying it 
to the remainder shell, we arrive at results for the first n 
eigenvectors in the N dimensional Hilbert space. 

It is obvious that this procedure introduces a new ele- 
ment of approximation, namely we partially neglect part 
the matrix elements mixing the higher shells (part of it 
is retained as the starting No basis vectors at step I will 
contain components from the previous I — 1 shells). In 
terms of the parameters there are two sources of error: 
setting i) k > 1 and ii) n ^ No- E.g. for n = No the 
procedure would be exact when k = 1, i.e. if the comple- 
tion is done in one step. The error resulting from these 
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sources is difficult to control. Nevertheless, our experi- 
ence shows that by choosing the starting dimension N n 
and the number of steps k well, the results are very much 
improved. For low values of the total dimension N, a di- 
rect comparison to the exact spectrum in the truncated 
space is possible and this was carried out as a preliminary 
study to assess the reliability of the method (see later) . 

In principle, once the cutoff is high enough one can 
apply a perturbative renormalization group to extrapo- 
late to infinite cutoff [32, 38-40]. This predicts that for 
a quantity M. with a scaling dimension y, the cutoff de- 
pendence would be 

M(e cut ) = Mac + A(e cut y y + . . . (4.1) 

Our numerical data shows that this scaling region has 
not yet been reached for the diagonal form factors using 
up to approximately 500000 states (e cut = 36). However, 
we observed that the more naive method of using the 
first order formula (4.1) and letting y to be a free fitting 
parameter works very well, and we always obtained an 
exponent y close to 2. 

This is explained by the following considerations. A 
simple guess for the next correction is 

M(e cut ) =Moo+ A(e cut )-y + B(e CMt )"^ (2 " 2A) + • • • 

2 2A (4 ' 2) 

since the coupling has a dimension of [energy] . In- 
deed, this is found to work surprisingly well. For the 
operator whose matrix elements are considered here one 
has y = 2 A, so that (4.2) reduces to 

M(e cut ) = Moo + A(e cut y 2A + B(e cut y 2 + ... 

and it turns out that in all of our calculations the cut- 
off dependence is dominated by the second exponent. In 
some cases, keeping the lower order term produces a false 
numerical fit, so sometimes its omission necessary in or- 
der to get a meaningful extrapolation. 

In Fig. 4.1 numerical evaluations of the vacuum ex- 
pectation value at £ = I are displayed. The NRG flow 
of |(0)| was initiated at e cu t ; o = 20 (5614 states). For 
e cut = 23 (16032 states) TCSA calculations with and 
without NRG were also carried out. The difference be- 
tween the two sets of data points is small. Note that 
using a step size of 100 states produced data differing 
from the displayed NRG-TCSA points only in the 6th 
digits. Results for e cu t = 31 (181046 states) and two- 
exponent extrapolation (4.2) using data up to e cut = 31 
are also shown. 

In conclusion, utilizing the NRG (and RG) improve- 
ment it is possible to obtain a satisfactory agreement 
between the theoretical and TCSA vacuum expectation 
values. Based on this preliminary study, it is also rea- 
sonable to believe that the NRG error remains small for 
other diagonal matrix elements, thus the NRG improved 
TCSA should be reliable. 
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Figure 4.1: TCSA evaluation of the vacuum expectation value 
at £ = 2/7, with and without NRG, compared to the exact 
result (2.3). An extrapolation using (4.2) is also shown. 

5. RESULTS 
5.1. Diagonal scattering 

Before presenting the more interesting results corrobo- 
rating our conjecture (3.2) we show some data confirming 
(3.1), which also serve to test the accuracy of the NRG- 
TCSA method. 

This case was already investigated in [7]. However, in 
that work we found that the agreement between theory 
and TCSA results was rather unsatisfactory, especially 
for £ = 2/7. In contrast, the NRG greatly improves the 
situation, as shown in Fig. 5.1 where two and three soli- 
ton diagonal matrix elements corresponding to the lowest 
lying states in the appropriate sectors are plotted. The 
data are parametrized by the compactification radius of 
the bosonic field defined by 

and related to £ as 



s 2R 2 - 1 

£ = 2/7 and £ = 50/311 correspond to R = 1.5 and 
1.9, respectively. For the two soliton data we employed a 
second order extrapolation (4.2) using NRG-TCSA data 
up to e cu t = 34 (corresponding to approx. 277 thousand 
states). For the three soliton data an extrapolation was 
used omitting the e~ 2 t A term and NRG-TCSA data up 
to e cu t = 37 (corresponding to approx. 140 thousand 
states). The agreement is quite good even in the less 
trivial three-soliton case, where six soliton form factors 
are encountered on the bootstrap side. 
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* ++<|l/2,-l/2)|0|{l/2,-l/2)> ++ (R=1.9) 

0.8- ■ ++ <{3/2,-l/2||0|{3/2,-l/2|> ++ (R=1.9) 

A ++ <|l/2,-l/2)|0||l/2,-l/2)> ++ (R=1.5) 




(b) 

Figure 5.1: Diagonal matrix elements with (a) 2 solitons and 
(b) 3 solitons, all of them of the same (+) charge. (Here and 
in the subsequent figures / denotes the matrix element values 
in units of M.) 



5.2. Non-diagonal scattering 

Now we turn to the case of non-diagonal scattering. In 
Fig. 5.2 a large number of comparative data is shown 
for diagonal matrix elements in the Q = sector and 
compactification radii R = 1.5, 1.7 and 1.9 (correspond- 
ing to £ = 2/7, 50/239 and 50/311, respectively). All 
the NRG-TCSA data are improved by second order RG 
extrapolation (4.2) using measurements up to e cu t = 31 
(approx. 180 thousand states). The agreement is again 
very convincing. Apart from the consistently small sys- 
tematic errors there are a few outlier points due to level 
crossings (for a discussion see [7]), and also the break- 
down of the extrapolation in larger volume is seen. The 
latter could again be remedied by omitting the e~^ t A term 
from (4.2). 



6. CONCLUSIONS 

In this paper we discussed diagonal matrix elements 
between multi-soliton states in sine-Gordon theory. Be- 
cause of the strong cutoff dependence of diagonal matrix 
elements observed in previous works [5, 7], we first had 
to improve the TCSA method by taking into account 



more states from the conformal Hilbert space. This was 
accomplished by implementing the numerical renormal- 
ization group method proposed by Konik and Adamov 
[32], and then improving the resulting NRG-TCSA data 
further by cutoff extrapolation, using the TCSA renor- 
malization group introduced in [32, 38]. 

The theoretical description for diagonal matrix ele- 
ments in states with diagonal scattering was known from 
[2], and so we could use expectation values computed in 
two-soliton states and three-soliton states, with all soli- 
tons having positive topological charge, to verify the ac- 
curacy of the method. In this way we demonstrated con- 
clusively that the discrepancies observed previously [7] 
were really due to cutoff effects. In fact, in our previous 
work we did not even attempt to show data for diago- 
nal matrix elements with Q = +3 three-soliton states 
because of the large inaccuracy, but the NRG-TCSA 
method allowed for a satisfactory comparison of these 
as well. 

The main theoretical result of this paper is a conjec- 
ture for diagonal matrix elements in states with non- 
diagonal scattering, formulated in Section 3. For the 
already known case of diagonal scattering, this conjec- 
ture reduces to the well established formula (3.1), first 
obtained in [2]. Therefore we verified the first nontriv- 
ial case, which was diagonal matrix elements in Q = 
two-soliton states, i.e. states containing a soliton and 
an antisoliton. The conjecture did fit the data in a very 
convincing way. 

It would be interesting to provide further verification 
with Q = +1 three-soliton states; we tried to identify 
such states in the spectrum, but due to the density of 
the spectrum and the exponential finite size corrections 
for small volume, this proved elusive so far. There are 
certain possibilities to improve this situation: using the 
exact NLIE description [41-45] of the finite volume lev- 
els instead of the approximation provided by the Bethe- 
Yang equations (2.13), and implementing further im- 
provements to the TCSA method. However, this is out- 
side the scope of this work, and we hope to return to this 
problem at some point in the future. 

Nevertheless, even by restricting ourselves to the re- 
sults for which we have sufficient numerical evidence, we 
have a complete description of all finite volume form fac- 
tors below the three-soliton threshold. This makes it pos- 
sible to evaluate the one-point and two-point functions 
including all corrections below the three-soliton thresh- 
old using the methods developed in [2] and in [18], re- 
spectively. We must stress that while here we considered 
sine-Gordon theory, the results are expected to be valid 
for general integrable models, such as the 0(3) nonlin- 
ear a model, and are therefore potentially applicable to a 
range of field theoretic and condensed matter problems. 
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R — 1.5, even states R — 1.5, odd states 




R — 1.7, even states R — 1.7, odd states 
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Figure 5.2: Diagonal matrix elements in even and odd soliton-antisoliton states 
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